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Abstract 

We discussed exact solutions of the Schrodinger equation for a two-dimensional parabolic confine- 
ment potential in a homogeneous external magnetic field. It turns out that the Jujo-electron system 
is exactly solvable in the sense, that the problem can be reduced to numerically solving one radial 
Schrodinger equation. For a denumerably infinite set of values of the effective oscillator frequency 
uj — \/Wq + (u> c /2) 2 (where ujq is the frequency of the harmonic confinement potential and lo c is the 
cyclotron frequency of the magnetic field) even analytical solutions can be given. Our solutions for 
three electrons are exact in the strong - and the weak correlation limit. For quantum dot lattices 
with Coulomb-correlations between the electrons in different dots exact solutions are given, provided 
the lattice constant is large compared with the dot diameters. We are investigating basic physical 
properties of these solutions like the formation and distortion of Wigner molecules, the dependence 
of the correlation strength from ujq and ui c , and we show that in general there is no exact Kohn-Sham 
system for the semi-relativistic Current-Density-Functional Theory. 

keywords: quantum dots, exact solutions of Schrodinger equation, Wigner molecules, Current-Density- 
Functional Theory 



1 



1 Introduction 



Exact (and somtimes analytical) solutions of the Schrodinger equation for realistic few- 
electron models of the quantum-dot-type provide a lot of unique opportunities. The 
physical essence of basic physical notions like the formation and distortion of Wigner 
molecules (WMs) and the consequences of inter-dot electron correlations in dot lattices 
can be understood more easily than with numerical brute-force approaches, which provide 
no formulas but only data. Moreover, these phenomena can be monitored over a wide 
range of external parameter values, which allows us to tune the system continuously 
between the weakly and strongly correlated regime. Intermediate and strong correlations 
are of particular interest, because for weak correlations there are a lot of mean-field 
approaches available. In particular we can easily see what the basic difference between 
strong correlated systems with low density (Wigner crystal) and high magnetic fields 
(fractional quantum Hall systems) is. Further, exact solutions allow us to check the 
precision of approximations like Hartree-Fock and Density Functional approaches and to 
reveal their weak and strong points. Unlike the comparison of approximate solutions 
for real systems with experiments, this approach has the advantage that all physical 
quantities (including those which are not experimentally accessible) can be considered, 
there are no experimental side effects, which obscure the comparison, and there are no 
discrepancies due to differences between the model and the real system. Last but not least, 
basic mathematical assumptions about the structure of solutions (like non-interacting v- 
representability in current density functional theory), which cannot be proven for general 
systems, can sometimes be rejected for special systems. 

This paper is not a comprehensive review on quantum dots and quantum dot molecules 
and lattices. In particular it does not describe the approaches and the results from nu- 
merical diagonalizations in a complete set of basis functions, quasi-classical approaches 
for the Wigner limit, as well as quantum Monte Carlo -, density matrix renormalization 
group -, current density functional -, and Hartree-Fock approaches. Each of these meth- 
ods warrants a separate review (see e.g. Ref.s [1-4] and references therein). Instead, this 
paper is focused on those systems which can be solved exactly or analytically, albeit the 
more complex sytems only in some limits for the external parameters. 

A problem with some exactly solvable models is that they have to be sufficiently simple 
and the question is whether all their features are shared by real systems of greater com- 
plexity. Therefore, all approaches complement one another and they should be pursued 
in parallel. Apart from this aspect, our models for the two and three-electron quantum 
dots and lattices from two-electron dots are already interesting on its own. 
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2 Specification of the model and exact solutions 



2.1 Model Hamiltonian 

We consider a two-dimensional (2D) two-electron system (with Coulomb interaction) in a 
harmonic scalar potential v ext (r) = (1/2) lo 1 r 2 and a magnetic field B = B e z represented 
by the vector potential (in symmetric gauge) A ext (r) = (1/2) B x r = (1/2) B r e a . We 
introduced cylinder coordinates (r, a, z) with the cylinder axis perpendicular to the plane, 
to which the electron motion is confined. The Hamiltonian reads 
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where H spin = g* s « ' B, and atomic units U — m — e — 1 are used throughout. This 
is a widely used effective Hamiltonian model for a two-electron quantum dot. 

2.2 Exact solutions of the Schrodinger equation 

The Schrodinger equation with the Hamiltonian ([1]) can be solved not only by reduction 
to the numerical solution of an (ordinary) radial Schrodinger equation [20], but even 
analytically for a discrete, but infinite set of effective frequencies uj = ^Ju 2 + (u c /2) 2 [5] , 
where we introduced the cyclotron frequency u> c = B/c. 

If we introduce relative and center of mass (cm.) coordinates 

r = r 2 -ri , R=^(r! + r 2 ) (2) 

the Hamiltonian (TjQ) decouples exactly. 

H = 2 H r + - H R + H spin (3) 

The Hamiltonian for the cm. motion agrees with the Hamiltonian of a non-interacting 
particle in a magnetic field 
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and only the relative Hamiltonian contains the electron-electron interaction 



H r = — 
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where we introduced rescaled parameters ojr = Ar = 2A(R), uj t = \u>o, 

A-r — 2-^-(^") 

(the indices V and J R' refer to the relative and cm. coordinate systems, respectively). 
The decoupling of H allows the ansatz 

$ = £(R)y(r)x(si,s 2 ) , ( 6 ) 
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where x( s i, s 2) is the singlet or triplet spin eigen-function. 

The eigen-functions of the cm. Hamiltonian (J4]) have the form 

U M (R) _ e iM * 

t~7^ RMiRh 



M = 0,±l,±2,... , (7) 

where the polar coordinates of the cm. vector are denoted by (R, A) and the radial 
functions Um(R) and Rm{R) can be found in standard textbooks. 
With the following ansatz for the relative motion 

p imoL ,. ( r \ 

^ = - 7 = ^ , m = 0,±l,±2,... (8) 
V27T r 1 ' 1 

the Schrodinger equation H r tp(r) = e r (p(r) gives rise to a radial Schrodinger equation for 

u(r) 

where the polar coordinates for the relative vector are denoted by (r, a), uo T = ^lj, e r = 
e r — |mu; c , and uj c = B/c. The solutions are subject to the normalization condition 

oo 

/ dr\u(r)\ 2 = 1. The Pauli principle demands that (because of the different particle 

o 

exchange symmetry of the spin eigen-functions) the relative angular momentum m has 
to be even and odd in the singlet and triplet state, respectively. There is no constraint 
for the cm. angular momentum M following from the Pauli principle. Because of the 
orthogonality of the coordinate transformation, the above described solutions are eigen- 
functions of the total orbital angular momentum with the eigenvalue Ml = M + m. 

FigJT] shows that the modulus of the orbital angular momentum of the ground state 
(GS) grows stepwise with increasing magnetic field. This implies that the spin state 
oscillates between singlet and triplet [21]. The Zeeman term and quenching of the singlet 
state for higher magnetic fields is not included in FigQ] The cm. excitations are not 
included as well, because they have no impact on the character of the ground state. 

In Figf5]the magnetic field and the e-e-interaction are successively added to the levels 
in the confinement only. We observe that the magnetic field removes the degeneracy with 
respect to the sign of Ml and breaks the symmetry with respect to up and downward 
directed fields. Angular momenta which are parallel to the field (positive) produce an 
magnetic moment which is anti-parallel to the field. They have a large interaction energy 
with the external field which shifts the levels upwards. In the opposite case the shift due to 
this contribution is downwards. Without e-e-interaction the GS has always Ml = 0. The 
shift due to by the e-e-interaction is positive definite and it decreases with increasing \Ml\- 
This can be explained with the radial equation which determines the contribution 
of the relative motion to the total energy. For large Ml = m and small r the last (e-e- 
interaction) term is dwarfed by the second (centrifugal potential) term. For large r the 
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Figure 1: Total energy for fixed confinement frequency u>o = 1 versus cyclotron frequency lo c (i.e. 
magnetic field). The cm. system is always in the ground state with M = 0. The relative angular 
momentum m is varied. The vertical lines show where the total orbital angular momentum Ml = M + m 
of the ground state changes. S and T indicates whether the ground state is singlet or triplet. Thick lines 
indicate states which can be NIVR. 




Figure 2: Energy levels for fixed external fields versus total orbital angular momentum. The level of the 
GS is indicated. We started with noninteracting electrons in the confinement only, added the magnetic 
field and the e-e-interaction. (In the energy unit uj we used lo c = 2.5 for all three cases.) 
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third term (effective confinement) is dominating in all cases. In a classical picture this 
means that two electrons rotating with a high angular momentum are separated by the 
centrifugal force, so that the job of the e-e-interaction is already largely done and the 
addition of the e-e-interaction does not change much. It is this Ml- dependence of the 
e-e-interaction shift which moves the GS to smaller Ml with increasing uj c . 

2.3 Analytical solutions 

In [5,6] it has been shown that the radial Schrodinger equation has simple analytical 
solutions for a discrete, but infinite set of effective oscillator frequencies uj, the pattern of 
which for \m\ = 1 can be seen in Figj3l The patterns for all \m\ look qualitatively similar. 
All solutions have the following form: 

u(r) = r H+l p ( r ) e -\** , (10) 

where p(r) is a finite polynomial of degree (n — 1). For 'non-soluble' systems, the poly- 
nomial has an infinite number of terms. 




500 



1 /co 

Figure 3: Reduced energies (energy over effective oscillator frequency) versus inverse effective oscillator 
frequency for relative angular momentum |m| = 1. The crosses indicate solvable states. The lines are 
just a guide for the eye and they connect states with the same node number k. n is the same for all 
horizontal rows of crosses with the same ordinate. 

In this review we only provide the results for the simplest analytical solutions. For 
n = 1 there is only a solution for infinite u>, which we call an asymptotic solution, because 
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it is exact for uj — > oo 

1 = , p(r) = l (11) 
uj 

This solution agrees not only with the Laughlin model wave function (WF), if the latter 
is applied to N = 2 and expressed in terms of the coordinates used here, but it is also 
the exact solution for non-interacting electrons (electrons without Coulomb interaction) 
in relative- and cm. coordinates. The corresponding WF has no node and is a ground 
state. 

For n = 2 there is one finite-field solution 

I = (2 W + 1) , P{r) = l + ¥¥ L-. , (12) 

which is a ground state as well. 

For n = 3 there is one asymptotic solution, which is a first excited state, 

1 r 2 

- = , p r = 1 - j—. (13) 

uj (H + i) v ' 

and one finite field solution, which is a ground state. 
1 

- = 2(4|m|+3) , 



UJ 



j. - 



p( r ) = 1 + T-, ; r + ; r— : ; r (14) 

KJ (2|m| + l) 2(2 \m\ + 1)(4 \m\ +3) v ' 

The exact energies e of the relative coordinate systems corresponding to these eigen- 
functions can be obtained from one compact formula. 

(e — -muj c ) = (\m\ + n) uj (15) 

If we compare this result with the spectrum of a single electron in a quantum dot (e — 
^muj c ) = (\m\ + 2k + 1) uj, where k is the node number or degree of excitation, then it 
becomes clear that in both cases (e — \muj c ) is a integer multiple of uj. We can also say 
that an analytical solution exists for those value of \m\ and u), for which one state of the 
interacting system is degenerate with one state of the non-interacting system. 

A great deal of approaches for the 2D electron gas in external fields is based on model 
wave functions (WF). For inspirations and checks in the limit N=2 it would be desirable 
to have analytical solutions of the two-electron problem for a wide range of external field 
values. One idea is to use the special analytical solutions ffTO]) with one of the exactly 
solvable polynomials given in (I11H14I) as discussed in Ref. [6]. If we use a special exact 
solution for a finite interval of external potential values we assume that the polynomial 
is independent of uj in this interval and depends only on \m\. Each of the choices (fni 
[141) provides a different approximation to the exact solution. Now we check, what the 
precision of these choices over a wide range of external fields is. 
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Figure 4: Projection of analytical model wave functions with polynomials of degree n onto the exact 
WF. The logarithm of the deviation of the squared matrix element from unity is on the ordinate. Full, 
dashed and dash-dotted curves belong to log(u>o)=-2, 0, and +2, respectively (see also first number in 
the legend). Thick, medium and thin lines belong to n=l,2, and 3, respectively. 



8 



FigJH compares the precision of model WFs with different n. The negative poles in the 
curves for log(uo) = +2 and n=2 and 3 indicate the vicinity of those lu, which provide 
exact solutions. It is seen that the solution for n > 1 is everywhere better than the 
solution for n — 1 (the latter is exact for infinite fields), and that all solutions become 
exact for uj — > oo. This means that the ansatz for the N-electron system proposed in [6] 
is definitely more precise than the Laughlin ansatz, if both are applied to the 2-electron 
system. We want to mention that the logarithm of the relative error in the energies 
calculated with the model WFs shows qualitatively the same behavior as the projection. 
It is only smaller in magnitude. 



2.4 Exact densities 

With (EJ) and (IE]), we obtain for the total density 

n(r) =2 I dr' |$(r 

the general expression 

1 r2n 



n(r) 



2tt 2 Jo 



da / dr' 



Rm(^\ r 2 + -r' 2 + rr'cosa^j 



2 r ]2 
u m (r') 



(16) 



(17) 



Because we are interested in the ground state only, we can safely use the cm. state for 
M = 0: Ro(R) = 2\/ti exp(—u)R 2 ) which allows to do one integration analytically leaving 
us with 



n{r) 



7T 



-2uj r 2 



dr' e - {uj/2)r ' 2 I Q (2urr')\u m (r') 



18) 



where I n (x) are the modified Bessel functions. 

The general expression for the paramagnetic current density 



J P (r) 



dr' $*(r, r') V$(r, r') - $(r, r') V<F(r, r') 



is somewhat complicated. Therefore, we give here only the formula for M = 



f(r) 



e a m — e 

7T 



-2w T l 



dr 



e Q r{r) 



(19) 



(20) 



As to be expected, the paramagnetic current density is proportional to the total angular 
momentum, points in azimuthal direction e a , and the scalar j p (r) depends only on the 
distance r from the center and not from the azimuthal angle . 

Although both formulas (1181) and (1201) rely on the functions u m (r), which are solutions 
of (Q, the analytical behavior for r — > can be expressed in terms of two positive definite 
integrals. 



A 
A 2 



dr e-^' 2 > 2 
dr r 2 e~^ l2)r 



u m (r 



u m {r) 



(21) 
(22) 
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After power series expansion of I n (x), we obtain 



n{r) 



f(r) 



-2uH 



7T 



4 ^ 2 -2wr 2 

m e r 

7T 



A + A 2 w 2 r 2 
A + 77 A 



2 „2 



(23) 
(24) 

(25) 



For the origin this means that n(0) = 4uA /tt is always finite and j p {0) = always 
vanishes. On the other hand, the derivative of the density at the origin 4^(0) = 



vanishes, but the derivative of the paramagnetic current density ^(0) = m(4uj 2 /tt)A is 
finite, unless m — 0. Besides, there is a relation which does not involve the radial WFs 
explicitly. 



dr 



(0) = m oj n(0) 



The exact vorticity, which has the form 7(1") = e z 7(7*) reads in this limit 



7(r J 



m2w 1 



~2 ^2 2 1 

or — r z + 
^0 



(26) 



(27) 



As will be seen in Sect. 4, the limit r 
non-interacting i> representability. 



is decisive for our proof of the violation of 



3 Formation of Wigner molecules and correlation strength 

A illustrative classical picture for a WM in an environment with rotational symmetry is 
a rotating and vibrating electron molecule. For a two-electron system this is a dumbbell- 
like object. We are going to show that this configuration is a manifestation of strong 
e-e-correlations and it is formed in the limit of small ujq or large u c . The issue is: why is 
small ouq equivalent to large uj c although the exact WF and consequently all distribution 
functions depend only on the effective confinement frequency uj = + (lu c /2) 2 where 
both ujq and u c have qualitatively the same influence. In particular, we will point out, 
how strong magnetic fields can cause strong correlations. This is not the same mechanism 
as for weak confinement (see also Ref. [8]). 

For illustration we use the density (I16p . which provides the distribution of the electrons 
in space, and the pair correlation function 

g(r)=<^\J2S(r t -r 1 -r)\^> (28) 

i<j 

which determines the distribution of the distance between two electrons. Both quantities 
depend only from the modulus: n(r) = n(r) and g(r) = g(r). In our system the latter is 
given by the radial part of the relative WF alone. 

9(r) = MrW= [ 4^ (29) 
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Figure 5: Density (upper) and pair correlation function (lower) for several correlation strength, loo, oj c , 
and the angular momentum of the ground state are given in the legend. The scaling parameters are: 
f=0.0901(A), 1.0192(B), 0.6345(C), 14.5022(D), and 13.6472(E) 
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In FigOboth quantities are shown for a few typical parameter sets. We have chosen scaled 
axes (tilted quantities) which allow us to show and compare different extreme cases in 
the same picture. On the abscissas the r-coordinate is scaled f = r/f with the average 
radius of the dot f = (1/2) / d 2 r r n(r) and the ordinates are scaled in such a way that 
the norms are conserved: / d 2 f n(f) = 2 and / d 2 f g(f) = 1. For weak correlations n(r) 
is peaked in the region of lowest potential energy (center of the dot) and g(r) is spread 
out over the whole range of non-vanishing density, allowing in particular small distances 
between the electrons. This regime is realized in curves A and B applying to strong 
(uq large) or medium confinement and small and medium uj c . Strong correlations are 
connected with sharply peaked densities at non-zero r confining the electrons on a ring 
with radius f. The pair correlation function is sharply peaked at a distance 2f which is 
the diameter of the ring. This means that the electrons are localized at a ring and have 
virtually antipodal positions [?, Taut-2einB]hat agrees with the above mentioned calssical 
picture of a WM. This minimizes the e-e-interaction energy in the limits allowed by the 
potential confinement energy without enhancing the kinetic energy. In terms of external 
parameters this can be realized in two scenarios. 

i) If uj c is small (or zero), the confinement has to be weak (ujq small). This can be 
concluded from the comparison of curve B and D. For weak confinement the state is 
spread out widely and the density is low. Low density implies the dominance of the the 
e-e-interaction over the kinetic energy (last term in radial Schrodinger equation (Q versus 
the first term) and a state which minimizes the e-e-interaction. The condensed matter 
analog to this state is the Wigner crystal. All in all, strong correlations in systems with 
low densities are produced by the dominance of the the e-e-interaction over the kinetic 
energy. 

ii) If uj c is medium or large (and the confinement not too strong), then the angular 
momentum of the GS is strongly negative (see curve C and E). Large (modulus of) angular 
momentum means strong centrifugal potential (second term in radial Schrodinger equation 
(Q) which drives the electrons away from the center and produces the ring structure in 
the density. At the same time the remaining e-e-term maximizes the e-e-distance within 
the limits set by the density. This is most pronounced in curves C and E. The case of 
small cjo and large uj c (which is not shown because of its numerical difficulties caused 
by of the extremely strong angular momentum) is even more strongly correlated. The 
condensed matter analog to this state is the fractional quantum Hall state. The state of 
curve C shows strong correlations as well, but does not have a large diameter (f = 0.6345) 
and consequently a low density. This proves that low density is not necessary for strong 
correlation, but a magnetic field can do the job alone. All in all, strong correlations in 
high magnetic fields are mainly produced by the high modulus of the angular momentum 
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of the ground state. 

From these considerations it follows that in our family of systems a suitable dimen- 
sionless quantitative definition of the correlation strength^ can be set up by the mean 
square radius of the pair correlation function r 2 = f cPrr 2 g(r) and its half width A 2 r = 
/ d 2 r (r — f) 2 g(r), where r — f d 2 rr g(r), according to 

This means, small half width and a peak at large r produce strong correlations. FigJH] 
shows this quantity as a function of the external field parameters. The steps are caused by 
a change in the angular momentum of the GS. It is obvious that for small uj c a weakening 
of the confinement increases s corr , but an increase of the magnetic field is much more 
effective, if it is connected to an increase of the modulus of the angular momentum of the 
GS, what happens in the region where the steps are found. 

Matulis and Peeters [8] investigated the same issue using qualitative asymptodic ex- 
pressions for the wave function instead of exact solutions. Nevertheless the resulting 
trends are the same, although their visual picture is different. 

4 Violation of non- interacting v-representability of the exact so- 
lutions of the Schrodinger equation 

The exact solutions of the special system considered in this review can be used to show 
that (unlike generally assumed) an exact Kohn-Sham (KS) system in the framework of 
semi-relativistic Current Density Functional Theory (CDFT) can exist only in special 
cases. 

In Density Functional Theory (DFT) it can be shown that for the GS the external 
potential is a functional of the density [14] (see also textbooks [15,16] with more modern 
approaches) 

v ext (v) £ $ 5 n(r) (31) 

v- 1 

which would imply non-interacting -u-representability (NIRV) or the existence of an exact 
Kohn-Sham system for the GSs, if the interacting and the non-interacting systems would 
have a common set of ground state densities. 

In the presence of a magnetic field and for (semi-relativistic) Current Density Func- 
tional Theory (CDFT), the generalization of V" 1 for the ground state still exists, but 
Vignale and Rasolt [10, 11] just presupposed the existence of the generalization of C [17] 

1 Unlike in common quantum chemistry language, the correlation strength defined here comprises all effects beyond the 
Hartree approximation, in particular it includes the effect of exchange. 



13 



iog(s corr ) 




o 
o 

w, 1 

O 



0.5 



correlation strength s 



corr 




— log(co )=-2 

— Iog(co )= 

— log(co )=+2 



-2 -1.5 



-0.5 



0.5 



1 1.5 



iog(co c ) 



Figure 6: Correlation strength versus external field parameters in log-scales for all axes. The steps, 
corresponding to parameters where the angular momentum of the GS changes (see FigUJ), are fully 
resolved only in the lower panel, which shows curves for 3 discrete loq. 
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implying that NIVR and the existence of a KS scheme has not been proven. Capelle and 
Vignale [17], on the other hand, have shown that there can be several external potentials 
yext provide the same WFs and densities 

Vr 4 (r) \ 

V|^( r ) *^A/"(r) (32) 

/ 

where V ext (r) and 7V(r) represent both external potentials (v ext (v) and A ext (r)) and both 
densities (n(r) and j p (r)), respectively. Hence, C cannot exist anymore as an unique 
mapping and the question of NIVR cannot be answered in this way. However, the exact 
solutions of the special system considered in this review can be used to show that an exact 
Kohn-Sham system or NIVR can exist only in the following special cases. 

All those states at non-zero B can be NIVR, which are continuously connected to the 
singlet and triplet ground states at B — (see also FigJT])- in more detail: 

If the GS is a singlet (total orbital angular momentum Ml is even) both densities can 
be NIVR if the vorticity 7(r) = V x (j p (r)/n(r)) of the exact solution vanishes. For 
Ml = this is trivially guaranteed because the paramagnetic current density vanishes. 
The vorticity based on the exact solutions for the higher Ml does not vanish, in particular 
for small r. In the limit r — > this can even be shown analytically. 

If the GS is a triplet (Ml is odd) and we assume circular symmetry for the KS system 
(the same symmetry as the real system) then only the exact states with |Ml| = 1 can be 
NIVR with KS states having angular momenta mi = and \m2\ = 1. 
Without specification of the symmetry of the KS system the condition for NIVR is that 
the small-r-exponents of the KS states are and 1. 

The proof of the statement for the singlet state is extremely simple and will be given 
here. The other proofs and more detailed information can be found in [9]. The question 
is if the (in this case doubly occupied) KS-WF </?(r) = R(r) e 1 ^ can be chosen in such a 
way that the density and the gauge invariant vorticity of the non-interacting KS system 
and the exact solution agree. 

n exac t(r) = riKsir) = 2 [R{r)} 2 (33) 

lexacti?) = 7A's( r ) = (34) 

Eq. (I3"3"l) defines the real part of the KS-WF. On the other hand, the vorticity of a two- 
electron singlet KS state vanishes exactly irrespective of the special form of -R(r) and C( r )- 
Therefore, equation (I3"4"|) can only be satisfied if the vorticity of the corresponding exact 
solution vanishes as well. FigJT] shows that this is not the case, in particular for small 
r the violation is massive. Eq. fT2Tl) provides ■jexacti®) = 2 M L oj which shows that the 
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Figure 7: The vorticities for luq — 1 and a few typical cyclotron frequencies uj c where the state with 
negative Ml is the ground state. The sign of jexactir) agrees with the sign of Ml- 

'degree of violation' grows with growing u. The exact vorticity vanishes only for the state 
with zero angular momentum, which is the GS for small magnetic fields (see also FigJT]). 

5 Distortion of the three-electron 
Wigner molecule 

As shown above, the two-electron system is the simplest system which exhibits the phe- 
nomenon of the formation of WMs in finite systems. For three electrons there is another 
effect, namely a Jahn- Teller-like distortion of the WM [7] shown schematically in FigJHl 
which can be investigated using the solutions described above. The point is that in the 
strong correlation limit the three-electron system can be decoupled into three indepen- 
dent pairs, the Schrodinger equation for which agrees (apart from a renormalisation of 
the interaction parameters) with the Schrodinger equation for the relative coordinates in 
the two electron system [6,7]. 

The Hamiltonian of the three-electron system reads 

H = lt\l fa + -A(r 4 )) 2 + k 2 rfl + £ ^— (35) 

where the the Zeeman term is disregarded because it has no influence on the spacial dis- 
tribution of the electrons, but shifts only the energies. We consider the unitary coordinate 
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Figure 8: Schematic picture of the distorted three-electron Wigner molecule with total orbital angular 
momentum Ml = —1 (left) and Ml = —2 (right). The thin lines depict the undistorted WM. 



transformation from the original position vectors to new ones 
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where a= 1/3—1/ \/3 and b = 


1/3 


■f I/a/3- The corresponding inverse 



(36) 



provides for the difference coordinates in the e-e-interaction terms 



V3 (x - x fe ) 



(37) 



where k) = (1, 2, 3) and cyclic permutations, and X = | Y^=i x i is the center of mass 
(cm.) in the new coordinates. It is a special feature of this transformation that the latter 
agrees with the cm. R in the original coordinates. This transformation provides the 
equivalent Hamiltonian 



3 r 
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-A(x, 



1 2 2 1 1 

2 y/3\xi- 



XI 
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Next we have to observe that in the strong correlation limit the uncertainty of the cm. X 
is small compared with the expectation value of the new coordinates Xj (see the appendix 
of Ref. [7]). Therefore X can be considered as a small perturbation and in zero order in 
X the Hamiltonian fl38l) decouples into a sum of independent quasi-particle Hamiltonians 



H (o) 



3 

i=i 



The Schrodinger equation for the quasi-particles 

h <a,(x) = e q <p q (x) 



(39) 



(40) 



is similar to the Schrodinger equation for the Hamiltonian in the relative coordinates (151) 
and therefore it can be solved exactly. In terms of these solutions, the total energy is a sum 



17 



of quasi-particle energies and the total orbital eigenfunction is a product of quasi-particle 
functions. 



91,92,93 



(xi,x 2 ,x 3 ) 



E, 



'91,92,93 



<P Q1 (xi) • ¥> 32 (x 2 ) ■ ^ 3 (x 3 ) 



93 



(41) 
(42) 



where comprises all quantum numbers. We can consider the quasi-particles as electron 
pairs, wherby their WFs </? g (xfc) describe the distance between two electrons (see Eq. 
(|37p ). The crucial point in explaining the distortion of the WM is the generalised Pauli 
exclusion principle for the quasi-particles states, i.e., the rules which determine the allowed 
combination of quantum numbers in Eqs. (1411) and (1421) in order to guarantee the anti- 
symmetry of the WF under electron transposition. These rules depend on the total 
spin (S) and orbital angular momentum (Ml) and they rule for some configurations the 
agreement of all three quantum numbers g.; out (see [6]). This means that the electron 
distances in all three electron pairs cannot agree giving rise to a distortion of the WM. 

The final result is the following [7]: 
In the ground state the electrons in an WM form an equilateral triangle (as might be 
expected from naive reasoning) only, if the state is a quartet (S = 3/2) and the orbital 
angular momentum is a magic quantum number (Ml = 3 m; m = integer). Otherwise 
the triangle in the ground state is isosceles. For Ml = (3m + 1) one of the sides is longer 
and for Ml = (3m — 1) one of the sides is shorter than the other two. 

6 Coulomb correlations between Quantum dots 

Other systems where the above described two-electron solutions play a crucial role are two- 
electron quantum dot molecules and quantum dot lattices, where the Coulomb correlation 
between electrons in differents dots is taken into account in the Van der Waals approxi- 
mation [22] . This means that the diameter of the dots must be small compared with the 
distance between the dots and that the overlap between WFs of different dots should be 
negligible. Then the Coulomb interaction between the electrons at r nk = R° + and 
r n i k i = R°, + u n i k i in different dots centered at R° and R° ; with (n ^ n') can be expanded 
in second order (dipole approximation) as 



wmiltoinanshere the dipole tensor T(R) = (1/ R 5 ) [ 3 R o R — R 2 I ] has been introduced. 

In the following the arrangement and number of dots is arbitrary, but for simpler 
notations we consider only dots with two electrons each. The bare confinement potential 
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|(R°-R°,) + (u„. fc -u n , fc ,)| 
+ -(u„fc — u.„/fc/) • T(R° — R n ,) • (u n fc — u n ik>) + 
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can vary from dot to dot. If we introduce for each dot a cm. coordinate R n = R° + U n 
with U n = (l/2)(u n i + u n2 ) and a relative coordinate r n = r n2 — r ni = u n2 — u n i, then 
the total Hamiltonian decouples 



H = H cm ({Rn}) + Y, Hrel,n (r„) 
n 

into a collective Hamiltonian 
H, 



2 ^ 2 
P n + - A(U„ 

c 



and a sum of individual intradot Hamiltonians 
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The force constant tensor C njn / of the collective Hamiltonian and the effective confinement 
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Figure 9: Magneto-phonon dispersion in a rectangular quantum dot lattice with lattice constants oi/aa = 
2 for the critical interaction strength and for the three magnetic fields given in the legend. The abszissa 
gives the magneto-phonon frequency u> in units of the bare confinement frequency too of the dots. The 
wave vector (qi, 52) shown on the ordinate varies along the edge of the irreducible Brillouin zone. 



The spectrum of the intradot excitations from (I45l) can be obtained with the methods 
for single quantum dots and its general features have been described in the previous 
sections. The Hamiltionian (l4"4l) describes magneto-phonon excitations. If the strength of 
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the interdot interaction reaches a critical value, a magneto-phonon mode can become soft 
indicating a lattice instability. Such a case is shown in Figj9] for a rectangular pariodic 
lattice. Solutions for a selection of dot dimers and periodic lattices are given and discussed 
in some detail in Ref. [22]. 
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